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Gamow- Teller strength distributions and /3-decay half-hves in neutron-deficient Kr and Sr isotopes 
are investigated within a deformed quasiparticle random phase approximation. The approach is 
based on a selfconsistent Skyrme Hartree-Fock mean field with pairing correlations and residual 
separable particle-hole and particle-particle forces. A simple two-level model is used to mix the 
nuclear shapes into the physical ground state. Good agreement with experiment is found with 
shape mixing coefficients which are consistent with those obtained phenomenologically from mixing 
of rotational bands. 
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I. INTRODUCTION 



Neutron-deficient isotopes in the A = 70—80 mass region are known to be interesting examples where the equilibrium 
shape of the nucleus is the result of a critical interplay between various nuclear structure effects [H] . The evolution 
T-H of the nuclear shape in isotopic chains is rather involved due to the existence of shell gaps at nucleon numbers 34, 
', 36, 38, and 40. Thus, adding or removing a few nucleons may lead to shape transitions in neighboring nuclei. Shape 

■ coexistence is also present in this mass region, where competing prolate, oblate, and spherical shapes are expected 
(-H [ in the same nucleus at close energies. Calculations of the equilibrium configurations in this mass region have been 

I performed within different approaches, such as the configuration-dependent shell-correction approach with deformed 
Woods-Saxon potentials j2i] , selfconsistent deformed Skyrme mean- field calculations [3i, |4|] , Hartree-Fock-Bogoliubov 
calculations with the Gogny force relativistic mean-field calculations 0, or the selfconsistent complex excited 
VAMPIR approach f7|. 

The first experimental evidence of shape coexistence in neutron-deficient Kr isotopes was reported in Ref. Q 
from irregularities observed at the bottom of the rotational bands. Since theii, a great deal of data on the low-lying 
" ^ I excitation spectrum, including transition probabilities, have become available [S. [lOl. [Til. [T^ . One of the experimental 
^ I i indications revealing the existence of shape coexistence in even-even nuclei is the observation of low-lying 0+ excited 

■ states. Each of the 0+ states is interpreted as a ground state correspondiiig to a different shape. Experimental 
. evidences for 0+ shape isomers in Kr isotopes have been reported in Refs. [1, [lO, [HI- These states are connected by 
' electric monopole (EO) transitions, whose strength is related to the change in the rms radius of the nucleus between 

initial and final states and carries information about the change in deformation and the overlap of the wave functions. 
The relationship between i?0 strength and shape mixing has been investigated in Refs. In these references, 

^\ [ shape mixing and shape coexistence were analyzed in terms of a two-level model, which was shown to be a simple but 

■ successful method to interpret the phenomenology. In this model, the strength is related to the amount of mixing of 
^ ' configurations with different deformations in the physical states. This simple two-level mixing model has been also 

successfully applied to understand the low-energy spectra of neutron-deficient Kr isotopes [12| . In these isotopes, 
regular bands are observed at high spins, which is a characteristic of well deformed shapes, but for the lowest states 
the regularity is lost, which is interpreted as evidence for a perturbation from other close states. The anomaly at low 
spin observed in the systematics of the moments of inertia (see Fig. 3.37 in Ref. [IB]) has also been interpreted as 
an evidence for mixed ground states in the lightest Sr isotopes. This mixing predicts that EO transition strengths 
should be observable in neutron-deficient Sr isotopes. An experiment similar to the study of the ^^Rb decay to ^"'Kr 
[l6j would be very helpful for the ''^Y decay to ''^Sr. 

The decay properties of nuclei in this mass region have been also investigated both theoretically [13, [H, [H, US] 
and experimentally [2l|, [l^, [IJ . Deformation has been identified in those works as a relevant issue to understand 
/3-decay properties. In particular, the nuclear deformation is crucial to perform reliable calculations of nuclear reaction 
rates and /3-decay half-lives in nuclei involved in the rapid-proton capture (rp) process of relevance in X-ray burst 
scenarios [2^ [1^ . 

In this work we investigate the decay properties of neutron-deficient Kr (^^'^''Kr) and Sr (^^'^'^Sr) isotopes within a 
deformed Hartree-Fock (HF) with Skyrme interactions and pairing correlations in BCS approximation. Residual spin- 
isospin interactions are also included and treated in quasiparticle random phase approximation (QRPA) [l8l.l27ll28l.l29j. 
In Sec. II a brief review of the theoretical formalism is presented. Sec. Ill contains the results obtained within this 
approach for the potential energy curves, electric monopole strengths, Gamow- Teller (GT) strength distributions, and 
/3-decay half-lives. The results are discussed in terms of shape mixing using a two-level model. Sec. IV summarizes 
the main conclusions. 
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II. THEORETICAL FORMALISM 

In this section we summarize briefly the theory involved in the microscopic calculations. More details can be found 
in Ref. [l^ . The method consists in a selfconsistent formalism based on a deformed Hartree-Fock mean field obtained 
with Skyrme interactions, including pairing correlations in the BCS approximation. The single-particle energies, wave 
functions, and occupation probabilities are generated from this mean field. Two Skyrme forces are considered in this 
paper to quantify the theoretical uncertainties caused by the use of different effective interactions. One is the force 
Sk3 [30l |. which is one of the simplest and oldest parametrizations, and has been successfully tested against many 
nuclear properties in spherical and deformed nuclei. The other force is SLy4 [3l| . which is an example of one of the 
most recent parametrizations including selected properties of unstable nuclei in the adjusting procedure. 

The solution of the HF equation is found using the formalism developed in Ref. [32i |. assuming time reversal and 
axial symmetry. The single-particle wave functions are expanded in terms of the eigenstates of an axially symmetric 
harmonic oscillator in cylindrical coordinates, using twelve major shells. The method also includes pairing between 
like nucleons in BCS approximation with fixed gap parameters for protons and neutrons, which are determined phe- 
nomenologically from the odd-even mass differences through a symmetric five term formula involving the experimental 
binding energies [s^ . 

The energy surfaces are analyzed as a function of the quadrupole deformation. For that purpose, constrained HF 
calculations are performed with a quadratic constraint [3j|. The HF energy is minimized under the constraint of 
keeping fixed the nuclear deformation. Calculations for GT strengths are performed subsequently for the equilibrium 
shapes of each nucleus, that is, for the solutions, in general deformed, for which minima are obtained in the energy 
surfaces. Since decays connecting different shapes are disfavored, similar shapes are assumed for the ground state of 
the parent nucleus and for all populated states in the daughter nucleus. 

To describe GT transitions, a spin-isospin residual interaction is added to the mean field. This interaction contains 
two parts, particle-hole (ph) and particle-particle (pp). The interaction in the ph channel is responsible for the position 
and structure of the GT resonance [18, 28] and it can be derived consistently from the same Skyrme interaction used to 
generate the mean field, through the second derivatives of the energy density functional with respect to the one-body 
densities. The ph residual interaction is finally written in a separable form by averaging the Landau-Migdal resulting 
force over the nuclear volume, as explained in Ref. [l8| . The pp part is a neutron-proton pairing force in the — 1^ 
coupling channel, which is also introduced as a separable force [29[. 

= -2«:^j'^5^(-i)^p+p_K, p^^Y. {^Y''k)^\ ^) ■ (2) 

K irv 

The coupling strengths used in this work are x^gt ~ ^■^'^ MeV and kFqj, = 0.03 MeV. 

The proton-neutron QRPA phonon operator for GT excitations in even-even nuclei is written as 

where (a) are quasiparticle creation (annihilation) operators, ujk are the QRPA excitation energies, and X^^, Y^J^ 
the forward and backward amplitudes, respectively. For even-even nuclei the allowed GT transition amplitudes in the 
intrinsic frame connecting the QRPA ground state |0) [F;^^ |0) = 0] to one-phonon states \u}k) \}^tiK 1^) ~ I'^-R')]' 
are given by 

{iUK\^Kt^\0)^TM^'' , (4) 

where 
with 



3 



v's are occupation amplitudes (w^ = 1 — w^) and spin matrix elements connecting neutron and proton states with 
spin operators 

E^'^ = (HaK|7r) . (7) 

The solutions of the QRPA equations are found by solving first a dispersion relation of fourth order in the excitation 
ener gies iv. The technical procedure to solve the QRPA equations is described in detail in Ref. [1^. The Ikeda sum 
rule [3^ is always fulfilled in these calculations. 

The GT strength B^{GT^) in the laboratory system for a transition /^^^(O+O) — > I fK f{\^ K) can be obtained as 

B^(Gr±) = [{^K=o |aot±| 0)' 5{uoK=o -uj) + 2 {uok=i O)' 5{uk=i - w)] , (8) 

in [g\/4:TT] units. To obtain this expression, the initial and final states in the laboratory frame have been expressed in 
terms of the intrinsic states using the Bohr-Mottelson factorization [SB^. The effect of angular momentimr projection 
is then, to a large extent, taken into account. 

The /3-decay half-life is obtained by summing up all the allowed transition probabilities weighted with some phase 
space factors up to states in the daughter nucleus with excitation energies lying below the corresponding Q-energy, 



where D = 6200 s and 0.77 is a standard quenching factor that takes into account in an effective way all the correlations 
[s?! which are not properly considered in the present approach. Since Fermi contributions are very small, the bare 
results can be recovered by scaling the results in this paper for B{GT) and T1/2 with the square of this quenching 
factor. 

In (3~^ / EC decay, the Fermi integral f{Z,uj) consists of two parts, positron emission and electron capture. In this 
work they are computed numerically for each value of the energy, as explained in Ref. 38] . The inclusion of the Fermi 
strength B^{F^) — [gy/i-K] < wji+jO >^ in the half-lives becomes important only for nuclei with Z > N. 

The calculation of these contributions in our case shows that they are negligible in the N ^ Z + 2 isotopes (^*Kr,^^Sr) 
and represents only a few percent correction in the case N — Z (^^Kr,^®Sr). 



III. RESULTS 



A. Potential energy curves and electric monopole strength 

The potential energy curves corresponding to the force Sk3 (SLy4) for the isotopes '^2,74-|^j. g^^^ 76,78gj. ^g^j^ seen 
in Fig. [I](l2]). We obtain an oblate ground state in ^^Kr with a prolate local minimum at higher energy. In the case of 
^^Kr, the oblate and prolate solutions are practically degenerate in the Sk3 case, while SLy4 predicts an oblate shape 
in the ground state. The Sr isotopes show for Sk3 a prolate ground state and a second oblate (spherical) minimum 
for ^^Sr (^^Sr). In the case of SLy4 the oblate and prolate solutions are practically degenerate in ^^Sr. The spherical 
shape becomes the ground state in the case of ^*Sr. These results are in qualitative agreement with similar results 
obtained in this mass region from different methods 0, d, [3j S B 0) • In all of these studies "^^Kr is found to be 
oblate in its ground state with a prolate solution at higher energy. '^''Kr is found to exhibit prolate and oblate solutions 
close in energy. Both ^^Sr and ^*Sr show prolate deformations in their ground states with oblate and spherical minima 
at higher energies, respectively. These features are also in agreement with experiment, where similar shape changes 
have been observed in this region [12, . 

As it was mentioned in the Introduction, the two-level mixing model has been successfully used to interpret the low- 
lying excitation spectra of neutron-deficient Kr isotopes measured in Coulomb excitation experiments [12l | . Coulomb 
excitation is a suitable method to distinguish between different shapes of the nucleus and to verify the shape coexistence 
scenario in light Kr isotopes. The analysis carried out in Ref. [T^ was based on the assumption of two regular 
rotational bands with collective intra-band transitions and zero matrix elements between the intrinsic states of the 
different bands. It confirms the prolate-oblate shape coexistence scenario in ^''Kr, with an extracted squared mixing 
amplitude for the prolate configuration A = 0.48. The squared amplitude is only A = 0.10 (mostly oblate) in ^^Kr, as 
extracted phenomenologically in Ref. [ll| . 

The electric monopole operator can be expressed in terms of single-nucleon degrees of freedom as 
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Tm = ^ 



ekrl . 



(10) 



The diagonal matrix elements of this operator give information about nuclear radii. The nondiagonal matrix elements 
give EO transition amplitudes. If we consider the mixing of two O"'' intrinsic deformed configurations, the mixed or 
physical states (ground state and excited state) can be written as a linear combination of them, 



0+,) ^ \/r^|o+,„i)-%/A|o+,) , 



(11) 



then, 
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(12) 



The cross term is neglected because the wave functions are mainly localized at different points in deformation space 
[T^ . This expression can be written in terms of the square of the difference between the rms radii of the states which 
are mixed. Using an expansion of the radii in the deformation parameter (3 up to third order, this expression can also 
be written as UM Eli 



A(l 



V/- pr 



/^obl) 



(13) 



As one can see, if there is no mixing, the -EO strength becomes zero and it is largest for maximal mixing. If the 
deformations P are known, one can extract from these expressions the mixing of the two deformed configurations. 
Dynamical effects beyond mean- field approach have been considered in various works 0, [13] ■ The results obtained 
from the Generator Coordinate Method and those from the phenomenological two-level model are compared in the 
mass region A — 190. The monopole strengths calculated with both methods agree quite well [4^ . 

Fig. [3] contains the electric monopole strength p^{EO) obtained with the force Sk3 from Eq. (I13p as a function of 
the mixing parameter A. The solid (dashed) lines correspond to the results obtained with (without) the terms in 
Eq. (|13D . terms that are neglected in many works. Horizontal lines in ''^Kr and ''^Kr are the experimental values 
from |11| and @, respectively. One can see that the experimental values in the Kr isotopes are compatible with 



A = 0.1 — 0.2, which is consistent with the mixing extracted in Ref. 
the mixing extracted in Ref. [l3 for ^^Kr (Aexp = 0.48). In the case of Sr isotopes 



11| for 
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Kr (Acxp = 0.1), but inconsistent with 
^(i?0) strengths are larger than 



in the case of Kr isotopes. Assuming that the experimental strengths are similar to those for Kr isotopes, this would 
indicate very weak mixing, which is consistent with the strong prolate component expected in these isotopes. 



B. Energy distributions of the Gamow- Teller strength 

In previous works [H, [2^, we have studied the sensitivity of the GT strength distributions to the various 
ingredients contributing to the QRPA-like calculations, namely to the NN effective force, to pairing correlations, 
to deformation, and to residual interactions. We found different sensitivities to them. In this work, all of these 
ingredients have been fixed to the most reasonable choices found previously. Here, we mainly discuss the mixing 
of different shapes needed to reproduce the available experimental information on the energy distributions of GT 
strengths. Experimental information on GT strength distributions is available for ^^Kr [2l|, ^"^Kr and ^^Sr [23| . 
Measurements on ^*Sr are being presently analyzed It should be mentioned that data for ^"^Kr, ^^Sr, and ^^Sr 
were taken using the total absorption gamma spectroscopy (TAgS) technique that avoids systematic uncertainties 
related to the so-called Pandemonium effect associated with the high resolution techniques, used for example in the 
case of ^^Kr [2l[. Thus, while the GT strength has been observed over most of the Qec window in the former cases. 
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in the latter case, only the GT strength below 2 MeV was extracted, which is still far from the Qec = 5.04 MeV 
energy. 

Figures [4] and [5] show the cumulative GT strength distributions as a function of the excitation energy in the daughter 
nucleus. Data are compared with Sk3 and SLy4 calculations corresponding to the various equilibrium configurations 
and to the adopted mixing. A quenching factor 0.77 is included in these calculations. When the ground state of the 
parent nucleus is given by a superposition of prolate and oblate shapes (Eq. (fTT|) ). the final 1+ states reached by the 
action of the GT operator are either jlp^.^j > or jl^j^j > at excitation energies given by Wproi and ojobi, respectively. 
Thus, 

B^GT) ^XJ2 (lproi> ^Proi \GT\ 0+„i)' <5K,oi - ^) + (1 - A) ^ (l+ c^obi \GT\ 0+ ^ J^bi - cj) . (14) 

^prol ^obl 

This is certainly a simple approach that gives us only a limited insight into the more involved problem of configuration 
mixing, but it is a first step in this direction toward the effect of mixing in the Gamow- Teller strength distributions. 

The mixing coefficients considered in these plots correspond to the experimental mixing obtained in Refs. [Ill) 
for the Kr isotopes. In the case of ^^Sr no mixing is plotted since the prolate configuration reproduces the data fairly 
well. In the case of ^*Sr, 50% mixing is plotted as a middle value useful to compare with future data. In the case 
of ^^Kr, the description of the GT strength distribution in the case of Sk3 is rather similar to both shapes. On the 
contrary, with SLy4, the oblate shape and in particular the mixing with a 10% prolate configuration describes the 
data fairly. Certainly, it would be interesting to extend the measurements up to the Qec window and compare the 
calculations in the whole range of energies. In the case of ^'*Kr, the experiment is better reproduced by the prolate 
shape with the two Skyrme forces. However, mixing with the oblate shape improves the results although 50% mixing 
seems to be very strong. 

The prolate configuration alone gives a good description of the GT strength distributions in ''^Sr with both Sk3 and 
SLy4 forces. The oblate shape generates a rather flat profile above 2 MeV that fails to account for the experimental 
profile of the strength distribution. The calculated strength distributions in ^*Sr show a very pronounced stepwise 
profile in the spherical case, as it corresponds to transitions between degenerate states. In the prolate case the strength 
is more fragmented. The total accumulated strength up to about 3 MeV is similar with the two forces, but in the 
case of Sk3, there is a strong transition just below Qec which doubles the strength. In the case of SLy4, this strong 
transition occurs a little bit above Qec = 3.76 MeV. 

These results also can be compared with other calculations such as those performed within a Tamm-Dancoff ap- 
proximation with Sk3 interaction p^ . In this reference, strengths contained in bins of 1 MeV were plotted and no 
further details within the Qec window were shown. Nevertheless, there is qualitative agreement with our results. 
More recently, results from the complex excited VAMPIR variational approach with Bonn potentials, on ^^Kr and 
^'^Kr, have been published [H, H^]. In the case of ^^Kr, the profile of the cumulative GT strength in Ref. [l^ presents 
a strong jump at around 1 MeV with a continuous increase elsewhere. This is at variance with ex peri ment which 
does not show this behavior. For '^^Kr, again a sudden increase of the strength at 1 MeV is found in [20] which does 
not show up in the experiment. The total strength is considerably underestimated by the calculations even though 
no quenching is considered in that work. 

C. Half-lives 

The half-lives calculated according to Eq.® can be seen in Table I. The results correspond to the different equi- 
librium deformations of the various isotopes, as well as to the mixed configurations obtained with similar mixing as 
those considered in Figs. S] and [5] for the GT strength distributions. Experimental Qec values have been used in 
these calculations. 

It should be stressed that the f{Z,uj) functions in Eq.® weight differently the strength depending on the energy. 
As a result, the half-lives are more sensitive to the strength located at certain excitation energies. One can see in 
Fig. [S]the f{Z,uj) functions corresponding to the /3+-decay, to the electron capture, and the total f{Z,uj) plotted 
versus the excitation energy of the daughter nucleus. The first thing to notice is that these factors are larger at lower 
excitation energies (higher energy of the /3-particle) and therefore, the contribution of the strength at low excitation 
energies to the half-lives is highly favored. One can also see that the /3+ component is the dominant contribution at 
lower excitation energy and it goes roughly like {Qp — E^x)^ ■ At higher energies, electron capture becomes dominant 
and it is indeed the only component between Qfj and Qec, where electron capture is allowed but positron emission 
is energetically forbidden. 

The results in Table I obtained for ^^Kr are once more compatible with predominantly oblate shape with a 10% 
prolate mixing [11|, although a little bit more mixing is favored (actually the half- life is reproduced with A = 0.3 in 
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the case of Sk3 and with A = 0.2 in the case of SLy4). In the case of ^^Kr one can see that the 50% maximal mixing 
considered also reproduces the experimental data reasonably well. In the case of the Sk3 force, A = 0.6 reproduces the 
half-life. The half-life of ^^Sr is compatible with a pure prolate shape, as in the case with the GT strength distribution. 
In the case of ^^Sr, a mixing of the spherical and prolate configurations in Sk3 reproduces the half-life, while in the 
case of SLy4 a prolate configuration is favored. 

IV. CONCLUSIONS 

In this paper /3-decay properties of neutron-deficient Kr and Sr isotopes are calculated within a deformed QRPA 
approach based on mean fields generated from selfconsistent Skyrme Hartree-Fock calculations. A simple two-state 
mixing model has been used to mix the intrinsic configurations into the physical states. It is shown that in some 
cases a single shape accounts for the main characteristics of the GT strength distributions and the half-lives. In other 
cases the data appear between the predictions of various shapes, demanding a more sophisticated treatment. Here 
we have considered a rough estimate of the mixing using the same coefficients extracted phenomenologically from 
Coulomb excitation experiments and have found that they reproduce the GT strength distributions and half-lives 
fairly well. Thus, we get a globally consistent picture when describing the isotope ^^Kr as a mainly oblate nucleus 
with small admixtures from a prolate configuration, the isotope ^^Kr as a more mixed nucleus with a dominant prolate 
configuration, the isotope ^^Sr as an almost prolate nucleus, and the isotope ''^Sr as a mixed spherical/prolate nucleus. 

The main objective in this paper has been to demonstrate that this approach is able to account for the main features 
of the decay properties of nuclei in this mass region, which are characterized by deformation as a key ingredient. This 
study provides additional and complementary indications in favor of shape coexistence in this mass region that is 
consistent with the information extracted from low-energy Coulomb excitation experiments. 

This approach cannot be pushed forward at a more quantitative description since the results are not only sensitive 
to the details of the present calculation but may also be sensitive to dynamical effects beyond mean-field approach 
not considered in this work. In this respect, it is worth pointing out that Multiple Coulomb Excitation data are 
beginning to yield sufficient numbers of E2 matrix elements for construction not only of quadrupole centroids but 
also fluctuation widths for the lowest few states in even-even nuclei Hi . 
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TABLE I: Half-lives {T1/2 [s]) corresponding to the forces Sk3 and SLy4 for the various shapes and prolate mixing A. In the 
case of '^*Sr the oblate results correspond actually to the spherical configuration. 



oblate prolate mixed (A ) exp. 


Sk3 












'■^Kr 


19.7 


14.0 18.9 (0.1) 


17.2 




^^Kr 


441.1 


1020.0 615.9 (0.5) 


690 




^«Sr 


4.2 


8.9 8.9 (1.0) 


8.9 




^«Sr 


93.8 


320.6 145.1 (0.5) 


150 


SLy4 


^^Kr 










18.2 


13.7 17.6 (0.1) 


17.2 




^*Kr 


718.9 


858.5 782.5 (0.5) 


690 




^«Sr 


5.2 


10.0 10.0 (1.0) 


8.9 






343.4 


165.9 223.7 (0.5) 


150 



Sk3 




-0.4 -0.2 0.2 0.4 -0.2 0.2 0.4 0.6 
P P 

FIG. 1: Potential energy curves for ^^'^*Kr (left panel) and '^®'^*Sr (right panel), obtained from constrained HF+BCS calcular 
tions with the Skyrme force Sk3. Beta is the quadrupole deformation. 




FIG. 2: Same as in Fig. [T]for SLy4. 
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0.15 




FIG. 3; Electric monopole strength p^{EO) with Sk3 as a function of the mixing parameter A, calculated from Eq. (|13|l with 
(solid lines) and without (dashed lines) (3'^ terms. Horizontal lines in ^^Kr and ^*Kr correspond to the experimental values 
from [ll|] and [3, respectively. 
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FIG. 4; (color online) Cumulative QRPA GT strength distributions with the force Sk3 as a function of the excitation energy 
in the daughter nucleus. The calculations correspond to the various equilibrium configurations and to the adopted mixing. 
Experimental data (black dots) are from [2l[ (^^Kr), [12] C^Ki), and Q C'^Sr). 
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SLy4 




FIG. 5: (color online) Same as in Fig. 14] but using SLy4 force. 



